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ABSTRACT 

We report the discovery of 3 new unbound hypervelocity stars (HVSs), stars traveling with such 
extreme velocities that dynamical ejection from a massive black hole (MBH) is their only suggested 
origin. We also detect a population of possibly bound HVSs. The significant asymmetry we observe 
in the velocity distribution - we find 26 stars with v^i > 275 km s~^ and 1 star with v^f < —275 
km s^^ - shows that the HVSs must be short-lived, probably 3-4 M0 main sequence stars. Any 
population of hypervelocity post-main sequence stars should contain stars falling back onto the Galaxy, 
contrary to the observations. The spatial distribution of HVSs also supports the main sequence 
interpretation: longer-lived 3 Mq HVSs fill our survey volume; shorter-lived 4 M0 HVSs are missing 
at faint magnitudes. We infer that there are 96 ± 10 HVSs of mass 3-4 Mq within R < 100 kpc, 
possibly enough HVSs to constrain ejection mechanisms and potential models. Depending on the 
mass function of HVSs, we predict that SEGUE may find up to 5 - 15 new HVSs. The travel times of 
our HVSs favor a continuous ejection process, although a ^120 Myr-old burst of HVSs is also allowed. 
Subject headings: Galaxy: halo — Galaxy: center — Galaxy: stellar content — Galaxy: kinematics 
and dynamics — stars: early-type 
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1. INTRODUCTION 

In 2005 we discovered the first HVS: a 3 Mq main 
sequence star traveling with a Galactic rest frame ve- 
locity of at least -1-709 ± 12 km s~^, many times the 
escape velocity of the Galaxy at its heliocentric dis - 
tance of 110 kpc (jBrown et al.ll2005l : lFuentes et al.ll2006D . 
This star cannot be explained by normal stellar interac- 
tions: the maximum ejection veloci ty from b inary dis- 
ruption mechanisms (jBlaauw 1961; Poveda et al. 196j) 
is limited to ~300 km s~^ for few M ^ stars (iLeonarc 



1991lll993H Tau ris fc Ta ken!jll998l:IPorte gies Zwartll200C 
Davies et a l. 2002; Gualandris ct al. 2005). Thus the ori- 



gin of the HVS must be tied to a more massive and com- 
pact ob ject. 

iHillg UMB) first predicted HVSs as an inevitable con- 
sequence of three-body interactions with a massive black 
hole (MBH). There is overwhelming eviden ce for a 4 x 10^ 
MfT) MBH at the cent er of our Galaxy (jSchodel et al.l 
[200l iGhez et al.l[2005h . Thus HVSs are probably stars 
ejected from the Galaxy by the MBH in the Galactic 
center. 

Furth er HVS disco v eries h ave followed the original dis- 
covery. iHirsch et al.l (|2005f ) reported a helium-rich sub- 
luminous O star leaving the Gala xy with a rest-frame ve - 
locity of at least -1-717 km s~^ lEdelmann et all (|2005l ) 
reported an 8 Mq main sequence star with a Galactic 
rest frame velocity of at least -1-548 km s~^, possibly 
ejected fr om the Large Magellanic Cloud. iBrown et al.l 
(I2006al lbl. I2007b[ ) designed a targeted HVS survey using 
the 6.5m MMT and Whipple 1.5m telescopes to measure 
the radial velocities of distant B-type stars. This strat- 
egy has worked remarkably well, yielding 7 HVSs and 
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evidence for a bound population of stars ejected by the 
same mechanism. Here we report the 3 newest HVS dis- 
coveries, and discuss additional evidence for the "bound" 
HVSs. 

The existence of HVSs has inspired broad theo- 
retical interest, and many testable predictions have 
emerged. It is clear that HVSs can be ejected by dif- 
ferent mechanisms: binary star encounters with a single 
MBH (pH^ fT988t TYu fc Tremai^ [2003t iBromlev et all 
|2006[) o r possibly with an interm ediate mass black hole 
(IMBH, I Gualandris fc Portegies" Zwart 2007), singl e star 
encounters with a binary M BH (Yu fc Trema inc 200^ 
ISesana et ail l2006t iMerrittI 2006; S esana et al., .2007ah 
or an in-sp i raling IMBH-MB H ( Gua landris et all 120051 : 
iLevinI 120061 : iBaumgardt et al.i,2006) . and single star en- 
counters with the cluster of stellar mas s black holes 
around the MBH (jO'Learv fc Loebl l2007f) . The dom- 
inant ejection mechanism remains unclear. Ejection 
rates depend on the number of stars scatt ered onto or- 
bits that intersect the MBH's " loss cone" (jPerets et al.l 
120071: iPerets fc Alexand"a1 120071 ). Interestingly, the dif- 
ferent ejection mechanisms result in different distribu- 
tions of HVS ej ection velocities an d eject ion rates (e.g. 
lYu fc Tremahd[2003l: ISesana et al.ll2007bD . 

HVSs probe a variety of characteristics of the Galaxy. 
The density of HVSs and their distribution of velocities 
tell us about the MBH's environment. The stellar types 
of HVSs t ell us about the types of stars orbiting near 
the MBH (IBrown et al.]|2006atlDemaraue fc Vira^ l2007t 
iKollmeier fc Gouldll2007l : iLu et al.ll2007D . The trajecto- 
ries of HVSs provide unique probes of the shape and ori- 
entat i on of the G alaxy's dark matter halo (|Gnedin et al.l 
120051 : lYu fc Madau 2007). A large sample of HVSs will 
be a new and powerful tool to investigate these issues. 
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Our HVS radial velocity survey is now 96.5% com- 
plete for faint B-type stars over 7300 deg^ of the 
Sloan Digital Sky Survey (SPSS) Data Release 5 (DR5, 
lAdelman-McCarthv et al.ll2007af ). Our survey provides 
strong evidence for a class of HVSs on bound orbits 
([Brown et alJ l2007bD matching exp ectations from theo- 
retical models (|Bromlev et al.ll20od: ISesana et aLll2007bD . 
The distribution of HVSs on the sky appears marginally 
anisotropic, as expected for a magnitude-li mited survey 
of HV Ss ejected from the Galactic center (jBrown et alJ 
l2007bD . 

Here we address the HVS's stellar type, their space 
density, and the history of HVS ejections from the Galac- 
tic center. In §2 we discuss the final target selection and 
spectroscopic identifications. In §3 we describe three new 
HVSs. In §4 we show that the observed HVSs must be 
short-lived objects, probably main sequence stars. In 
§5 we calculate the space density of HVSs, and predict 
HVSs discovery rates for some future surveys. In §6 we 
calculate HVS travel times, and discuss their ejection 
history. We conclude in §7. The new observations are in 
Appendix A. 

2. DATA 
2.1. Target Selection 

iBrown et al.l (|2006bl hereafter Paper I) describe the 
target selection for our survey of faint 17 < < 19.5 
HVS candidates. Briefly, we use SDSS photometry to se- 
lect HVS candidates with the colors of late B-type stars. 
B-type stars have lifetimes consistent with travel times 
from the Galactic center but are not a normally expected 
Galactic halo population. Our color selection introduces 
no kinematic bias and allows for stars at any velocity. 

In this paper we make two changes to the Paper I tar- 
get selection: 1) we select targets from the larger SDSS 
DR5 catalog, and 2) we include stars in a supplementary 
color-color region defined by 0.55 < {u' — g')o < 0.9 and 
-0.28 < (ff'-r')o < -0.25 (see Figured]). Our fuU target 
list contains 621 HVS candidates over 7300 deg^ of sky. 
We have observed 575 of these targets; we are thus 93% 
complete for the faint HVS candidates. The remaining 
4 6 targets are confined to the region 13^^25 < RA < 14^*5. 

iBrown et al.l (|2007bl hereafter Paper II) describe the 
target selection for our complementary survey of bright 
15 < g'o < 17 HVS candidates. We now select targets 
from the SDSS DR5 catalog. We identify a total of 892 
objects with B-type colors; 691 satisfy our Galactic lon- 
gitude and latitude cut. We have observed all 691 objects 
and are thus 100% complete for bright HVS candidates. 
These objects are spread over the identical 7300 dcg^ of 
sky as the faint HVS candidates. 

2.2. Spectroscopic Observations 

New observations of faint HVS candidates were ob- 
tained at the 6.5m MMT telescope with the MMT Blue 
Channel spectrograph on the nights of 2006 December 
21-26, 2007 March 18-20, and 2007 May 14-18. Wc op- 
erated the spectrograph with the 832 line mm^^ grating 
in second order and a 1.25" slit. These settings provide 
wavelength coverage 3650 A to 4500 A and a spectral 
resolution of 1.2 A. Exposure times ranged from 5 to 
30 minutes and yielded signal-to-noise S/N = 15 in the 
continuum at 4000 A. We obtained comparison lamp ex- 
posures after every exposure. 
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Fig. 1. — Color-color diagram showing the spectroscopic iden- 
tifications of the 1005 late-B spectral type stars (squares) and the 
244 white dwarfs (x's) in our survey. Target selection regions are 
shown for the MMT [short dashed line) and FAST {long dashed 
line) samples . If th e HVSs (stars) are main sequence stars, then 
IGirardi et al.l I I2004I ) stellar evolution tracks suggest they are 3 Mq 
(lower solid line) or 4 Mq stars (upper solid line). The average 
color uncertainties are illustrated by the errorbar on the lower left. 
Also marked are the 19 possibly bound HVSs in the MMT sample 
(triangles) and FAST sample (circles), and HVSl (plus sign). 



We also obtained new observations of bright HVS can- 
didates at the Whipple Observatory 1 .5m Tillinghast 
telesc ope with the FAST spectrograph (|Fabricant et al.1 
[1991 over the course of 11 nights spread between 2006 
September and 2007 May. We operated the spectrograph 
with a 600 line mm~^ grating and a 2" slit. These set- 
tings provide wavelength coverage 3500 A to 5400 A and 
a spectral resolution of 2.3 A. Like the MMT observa- 
tions, exposure times ranged from 5 to 30 minutes and 
yield S/N = 15 in the continuum at 4000 A. We obtained 
comparison lamp exposures after every exposure. 

We extracted the spectra using IRAF^ in the stan- 
dard way and measured radial velocities using th e 
cross-correla tion package RVSAO ([Kurtz fc Minkll 1998( 1. 
IBrown et al.l ([MOSiJ) describe the cross-correlation tem- 
plates we use. The radial velocity accuracy is ±11 km 
s^^ for the B-type stars. 

2.3. Spectroscopic Identifications 

Figure [T] plots the colors and spectroscopic identifica- 
tions of the 1,266 HVS candidates in our survey. We 
find that 1005 (79.4%) are stars of B spectral type and 
244 (19.3%) are DA white dwarfs. Curiously, the white 
dwarfs in our survey have very unusual colors for DA 
white dwarfs; these colors imply very low surface grav- 
ity. Indeed, a detailed analysis of the spectra reveal that 
one of the white dwarfs has a mass of 0.17 Mq, the lowest 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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mass white dwarf ever found (|Kilic et al.ll2007allB . Other 
remarkable o bjects include B sup ergiants in the Leo A 
dwarf galaxy (|Brown et al.ll2007al ) and one compact, ex- 
tremely metal poor galaxy with properties similar to the 
host galaxies of nearby ga mma ray bursts (|Kewlev et al.l 
IMFtI: (Brown et al.l[20073) . 

Our new observations include a B supergiant in 
the Sextans B dwarf galaxy. The star SDSS 
J095951. 18+05212 4.52 is locate d within the galaxy's 3' 
Holmberg radius (iMateol |1998) and has a +295 ±11 
km s^^ heliocentric radial velocity, consistent with the 
+295 ± 10 km s~^ systemic velocity of Sextans B 
(jFisher fc Tullvl Il975f ). The observed spectrum is that 
of a luminosity class I B supergiant. The (m — M) = 
25.6 distance modulus of Sextans B ijPiotto et al.lll99"^ 
iMendez et~an [2002) implies the star has My = -6.25, 
also consistent with a B I supergiant. We find no other 
B supergiants in our survey, suggesting that, to the depth 
of our survey, the census of star-forming dwarf galaxies 
in the Local Group may be complete over the region of 
sky covered by SDSS DR5. 

In Figure [T| we also plot iGirardi et all (|2002l |2004 
main sequence track for a solar abundance 3 star 
(the lower solid line) and a 4 Mq star (the upper solid 
line) . If main sequence stars fall in our survey, our color 
selection targets primarily 3-4 M0 stars. All the HVSs 
overlap the main sequence tracks within their l-a uncer- 
tainties {(7u'-g' = ±0.042, cTg'-r' = ±0.033). 

2.4. Radial Velocity Distribution and Possibly Bound 

HVSs 

Figure [5] plots the distribution of line-of-sight veloc- 
ities, corrected to the Galactic rest-frame Vrf, for our 
entire sample of 1005 B-type stars. We calculate mini- 
mum Galactic rest-frame velocities by removing the ra- 
dial component of the Sun's 220 km s~^ orbital motion 
and the radial component of the Sun's motion with re- 
spect to the Local Standard of Rest from the observed 
radial velocities (see Paper I). 

The tail of positive velocity outliers in Figure [2] is truly 
striking: we observe 26 stars with v^i > 275 km s~^ and 
only 1 star with v^i < —275 km s~^. We estimate the 
significance of the asymmetric velocity distribution by 
fitting a Gaussian distribution to the data. We itera- 
tively clip 3-cr outliers and find a +13 ± 3 km s"-'^ mean 
and a 106 ± 4 km s^^ dispersion. The lower panel of 
Figure [2] plots the residuals of the observations from this 
Gaussian distribution, normalized by the value of the 
Gaussian. Stars with velocities |wrf| < 275 km s^^ show 
low-significance deviations from a Gaussian distribution 
and are likely a normal halo population (Paper I). Inte- 
grating the tail of the Gaussian suggests we should expect 
7 stars with v^i > 275 km s~^. Yet we observe 26 such 
objects, including the unbound HVSs. There is less than 
a 10~^ probability of drawing 26 stars with Vi-{ > 275 
km from a Gaussian distribution with the observed 
parameters. Thus the observed asymmetry appears sig- 
nificant at the 5-(T level. 

We now ask whether the non-zero mean of the distri- 
bution is consistent with models of HVS ejections. If 
the 19 objects in excess of the Gaussian distr i bution 
Vrf > 275 km s"^ are HVSs, the iBromlev et all (|2006l ) 
ejection model for 3 Mq stars suggests that an additional 
11 HVSs have v^-f < 275 km s~^ in our survey. We ran- 
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Fig. 2. — Galactic rest-frame velocity histogram for all 1005 
B-type stars in our sample (upper panel). The best-fit Gaussian 
(thin line) has dispersion 106±4 km . The lower panel plots the 
residuals of the observations from the best-fit Gaussian, normalized 
by the value of the Gaussian. The strong asymmetry of stars with 
\v^f \ > 275 km is significant at the ~5(T level, and shows that 
the bound stars with +300 km s~^ are short-lived; we observe only 
1 star falling back onto the Galaxy near —300 km s~^. 



domly draw 30 objects from the IBromlev et al.l (j2006f ) 
HVS ejection model and add them to 975 objects drawn 
from a Gaussian distribution with dispersion 106 km s^^ 
and zero mean. The expected 30 HVSs shift the mean of 
the distribution by +2 ± 4 km s~^, consistent with but 
not entirely sufficient to explain the observed shift in the 
mean. 

In Paper II, we argue that the significant excess of 
stars with velocities around Vj-i ~ +300 km demon- 
strates a population of HVSs ejected onto bound tra- 
jectories. HVS ejection mechanisms naturally produce 
broad spe ctrum of ejection v elocities that include bound 
orbits (Br omlev et al.l [20061 ). We observe 19 possibly 
bound HVSs with 275 < v^t < 400 km s~^ in our com- 
pleted survey. The colors of the bound HVSs suggest 
they are a mixed population: some of the bound HVSs 
are significantly redder or bluer than the main sequence 
tracks (see Figure [T]). However, it is possible to pick 
a sample of ^12 bound HVSs, the number in excess of 
the Gaussian distribution, that overlap the color distri- 
bution of the unbound HVSs. Follow-up spectroscopy 
is therefore required to constrain the stellar type of the 
possibly bound HVSs. Because we do not know which of 
the possibly bound HVSs are true HVSs, we focus our 
attention on the well-defined sample of unbound HVSs 
in this paper. 

3. HYPERVELOCITY STARS 

3.1. Three New HVSs 

We discover 3 HVSs in our new observations. The 
objects are SDSS J094214.04+200322.07 (HVSS), 
SDSS J102137.08-005234.77 (HVS9), and SDSS 
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TABLE 1 
HYPERVELOCITY STARS 



ID 


Type 


Mv 


V 


Rgc 


I 


b 






tac 


Catalog 






mag 


mag 


kpc 


deg 


deg 


km s^^ 


km s^^ 


Myr 




hvsl 


B 


-0.3 


19.84 


111 


227.33 


-1-31.33 


840 


696 


145 


SDSS J090744.99+024506.91 


hvs2 


sdO 


+2.6 


19.05 


26 


175.99 


-1-47.05 


708 


717 


32 


US 708^ 


hvs3 


B 


-2.7 


16.20 


62 


263.04 


-40.91 


723 


548 


100 


HE 0437-5439^ 


hvs4 


B 


-0.9 


18.50 


82 


194.76 


-1-42.56 


611 


567 


125 


SDSS J091301.01+305119.8* 


hvs5 


B 


-0.3 


17.70 


45 


146.23 


+38.70 


551 


647 


60 


SDSS J091759.48+672238.3* 


hvs6 


B 


-0.3 


19.11 


78 


243.12 


+59.56 


626 


528 


130 


SDSS JIIO557.45+O93439.55 


hvs7 


B 


-0.9 


17.80 


56 


263.83 


+57.95 


534 


421 


110 


SDSS J113312.12+010824.9'5 


hvs8 


B 


-0.3 


18.09 


53 


211.70 


+46.33 


511 


429 


100 


SDSS J094214.04+200322.1 


hvs9 


B 


-0.3 


18.76 


68 


244.63 


+44.38 


628 


485 


120 


SDSS J102137.08-005234.8 


hvslO 


B 


-0.3 


19.36 


87 


249.93 


+75.72 


478 


432 


165 


SDSS J120337.85+180250.4 
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Il2006bl) 



J120337.86+180250.35 (HVSIO). HVS8 has a B9 spec- 
tral type, a -|-512±10 km s~^ heliocentric radial velocity, 
and a minimum velocity of +430 km s^^ in the Galactic 
rest frame. A sol ar metallicity 3 M(7) m ain sequence star 
has My ^ -0.3 (jSchaller et al.lll992[ ). This luminosity 
places HVS8 at a Galactocentric distance R — 53 kpc. 
The mass of the Galaxy within 50 kpc is 5.5 x 1 0^^ M© 
fWilkinso n'fc EvanslfToMISakamoto et al.ll2003f ). which 
implies that the the escape velocity at 50 kpc is ^^300 
km s"^ We conclude that HVS8 is very likely unbound 
to the Galaxy. 

HVS9 has a B9 spectral type, a -|-632±11 km s~^ helio- 
centric radial velocity, and a minimum velocity of -1-489 
km in the Galactic rest frame. Unlike the other HVSs 
in our survey, which are located 10° - 20° away on the 
sky from the nearest Local Group galaxy, HVS9 is 2?3 
from the Sextans dwarf galaxy. Any physical associa- 
tion with Sextans is very unlikely, howev er. Sextans is 
1320 ±40 kpc distant (Dolp hin et'an i2003') and has a he- 
hoccntric velocity of 224 ±2 km s'^ (Young 2000). Thus 
HVS9 is moving toward the dwarf galaxy with a relative 
velocity of 408 km s~^. Interestingly, HVS9 also has the 
reddest {g' — r') color of the HVSs; it is arguably the 
best blue horizontal branch (BHB) candidate among our 
HVSs. The equivalent width of its Ca 1 1 K li ne suggests 
low metallicity and the IClewlev et al.l ()2004f ) line-shape 
technique suggests low surface gravity. However, both 
measurements are notoriously uncertain at the hot effec- 
tive temperature of HVS9, (B - y) ~ -0.01. Whether 
HVS9 is a BHB star located at i? = 35 kpc or a 3 Mq 
main sequence star at i? = 68 kpc, its large velocity 
means it is certainly unbound. 

HVSIO has a B9 spectral type, a -f478 ± 10 km s^^ 
heliocentric radial velocity, and a minimum velocity of 
+432 km s~^ in the Galactic rest frame. HVSIO is the 
faintest HVS discovered by this survey, g' — 19.295 ± 
0.024, and is thus the most distant HVS known except 
for HVSl. Assuming it is a 3 main sequence star, 
HVSIO is ~84 kpc above the Galactic plane at its latitude 
5 = +76°. 

Table [1] summarizes the HVS properties. Columns in- 
clude HVS number, stellar type, inferred absolute mag- 
nitude My, apparent V magnitude derived from SDSS 
photometry, Galactocentric distance i?. Galactic coor- 



dinates {l,h), heliocentric radial velocity Uq, minimum 
Galactic rest-frame velocity Wrf (not a full space velocity), 
travel time estimate from the Galactic Center taCi and 
catalog identification. Velocities have changed slightly 
from previous work because we have obtained multiple 
observations of all the HVSs; here we report the weighted 
average of the velocity measurements. We include HVSl, 
HVS2, and HVS3 in Table [T] but exclude them from 
our analyses; the first three HVSs do not fall within the 
color/magnitude bounds of this targeted survey. 

We do not report errors in Table [T] because formal un- 
certainties are misleadingly small compared to the (un- 
known) systematic errors. For example, heliocentric ra- 
dial velocities are accurate to ±11 km s~^, but we have 
no constraint on the proper motion component of the 
rest frame velocity Wj-f. The luminosity estimates de- 
scribed below are precise at the 10% level for main se- 
quence stars, however the luminosity estimates could be 
over-estimated by an order of magnitude for post-main 
sequence stars. 

3.2. Constraints on HVS Binaries 

Approximately 25% of late B-type stars in the solar 
neighborhood are in binaries with periods < 10 days 
and with mass ratios greater than TO2/TO1 > 0.1 (jWolfil 
119781 ). Thus, in principle, 2±1 of the 8 late B-type HVSs 
(we include HVSl) could be spectroscopic binaries. Al- 
though HVSs are ejected as s ingle stars by standard HVS 
mechanisms. iLu et al.l (|2007D argue that the MBH binary 
mechanism can eject a compact HVS binary. Detection 
of a single HVS binary might thus provide compelling ev- 
idence for a binary MBH in the Galactic center (|Lu et aTl 
I200I . 

We have two spectroscopic observations of each of the 
7 HVSs in this survey, typically obtained a few months 
apart. In all cases, the radial velocities are identical 
within the accuracy of the measurements, ±11 km s~^. 
Thus it appears unlikely that the HVSs are compact bi- 
naries. We test this claim by comparing the observations 
with velocities drawn 10^ times from a typical 50 km 
semi-amplitude binary with random orbital phase and 
inclination. A Kolmogorov-Smirnov (KS) test finds a 
0.04 likelihood of drawing the observations from a bi- 
nary system with 50 km s~^ semi-amplitude. Additional 
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Fig. 3. — Summed HVS spectrum, created from the weighted 
average of our observations of HVSl and HVS4 - HVSIO, shifted 
to rest frame. The spectral type is that of a B8 - B9 star. The 
wavelength difference between the pair of Ca II K lines, one from 
the HVSs and one from the local interstellar medium that appears 
in the spectra of HVSS and HVS9, visibly indicates the large space 
motion of the HVSs, AA/A ~ 550 km s'^ . 



observations of higher accuracy are necessary to rule out 
radial velocity variations completely. 

We have also obtained three observations of HVSl 
spread over 3.5 years. The heliocentric radial velocities 
are 853 ± 12, 816 ± 14, and 832 ± 13 km s'^. Thus 
HVSl may have velocity va riations. However, H VSl is a 
slowly pulsating B variable (jFuentes et al.|[2006l ). Slowly 
pulsating B vari ables exhibit radial velocity ampHtude s 
of -20 km s-i (jAerts et al.lll999HMathias et al.ll2001f ). 
consistent with the observations. Although our current 
observations are not conclusive, it is unlikely that any of 
the known HVSs are compact binary systems. 

3.3. HVS Template Spectrum 

Summing our HVS spectra allows construction of a 
high S/N template of the HVSs. Figure [3] displays the 
weighted average of our observations of HVSl and HVS4 
- HVSIO, shifted to the rest frame. The total integration 
time is equivalent to 4 hours on the 6.5m MMT tele- 
scope. The strongest spectral features in Figure [3] are 
the hydrogen Balmer lines, visible from H7 to H17. 

It is interesting to discuss spectral classification in 
light of the weak spectral lines visible in Figure \3\ The 
strengths of Mg 11 4481 and the Si 11 4130 blend rel- 
ative to He I 4471 and He i 4144, respectively, indicate 
that the summed HVS spectrum must be later than a B7 
spectral type. Notably, Ca 11 3933 and He i 4026 have 
similar equivalent widths, indicating that the summed 
HVS spectrum has a spectral type of B8 - B9. This late 
B spectral type is consistent with our color selection and 
our previous spectral classifications. 

Amusingly, we also see two Ca 11 K lines in Figure [3l 
One Ca 11 K line is from the atmospheres of the HVSs 
at 3933 A, and the other is from the local interstellar 
medium in the direction of HVS8 and HVS9. The in- 
terstellar Ca II K line appears shifted to 3926 A in the 
rest- frame of HVS8 and HVS9. The difference between 
the two Ca 11 K lines, AA/A — 550 km s^^, visibly indi- 
cates the large space motion of the HVSs. 



4. NATURE OF THE HYPERVELOCITY STARS 

We use the observed velocity distribution (Figure [5]) to 
show that HVSs must be short-lived. Using stellar evo- 
lution arguments and orbit calculations, we argue that 
the HVSs are likely main sequence stars. 

4.1. The Asymmetric Velocity Distribution 

The HVSs are moving away from us with large positive 
radial velocities 400 < v^i < 700 km s^^. As stated 
earlier, the escape velocity of the Galaxy at 50 kpc is 
~300 km s-\ or 360 km s^^ in the lBromlev et al.l ()2006D 
model (see Figure H}. Therefore the absence of objects 
with < —400 km s~^ is no surprise: unbound HVSs 
will never come back. 

The significant excess of -1-300 km s~^ stars compared 
to —300 km stars shows that the possibly bound 
HVSs must be short-lived stars: every star traveling at 
-1-300 km s^^ must eventually fall back onto the Galaxy 
at rou ghly -300 km s'^. [Koilmcicr & Gould (2007) and 
lYu fc Madaui (|2007[ ) also make this point. The return 
time depends on the apocenter of the orbits, and ranges 
from 800 Myr to 2400 Myr for the range of distances and 
velocities we observe. Although a HVS's radial orbit may 
be altered by the a sphericity of the Galactic potential 
(|Yu fc Madaull2007l ). the effect on radial velocity should 
be minor at the distances probed by our survey. 

4.2. HVSs are Main Sequence Stars 

The stars we observe have the spectral types of B stars, 
but they may be main sequence or post-main sequence 
stars. After a main sequence star exhausts its hydrogen 
fuel, it evolves up the red giant branch, ignites helium, 
and becomes a blue horizontal branch (BHB) star. All 
stars with masses capable of igniting helium burning, in- 
cluding tur n-off stars in th e halo, spend of order 10* yr as 
BHB stars (|Yi et al.lll997l) . The ambiguity in identifying 
BHB stars arises because the horizontal branch overlaps 
the main sequence at the effective temperatures of our 
late B-type stars, T^s 10'' K. Thus we cannot discrim- 
inate BHB and main sequence stars by surface gravity 
alone. 

In Paper I, we use kinematics and metallicities to show 
that most of the B-type stars are consistent with a halo 
population of post-main sequence stars or blue stragglers 
in the halo. The HVSs are a different class of objects. 
Thus the question stands: are the HVSs 3-4 M0 main 
sequence stars or BHB stars? 

The colors we observe preferentially sample hot BHB 
stars with small hydrogen envelopes. Hot BHB stars 
are p resent in only 25% of globular clusters ()Lee et al.l 
|2007| ). One may wonder whether hot BHB stars 
are common in a Galactic cente r population with so- 
lar abundance (ICarr et all l2000t iRamirez et al.l l2000l : 
ICunha et al.l 120071 ). Theoretical evoluti onary tracks a l- 
low for BHB stars at all metallicities ()Yi et al.lll997f l. 
Hot BHB stars may, in fact, be more common at high 
metallicities because increased opacity may cause in- 
creas ed mass loss in the red giant branch phase (jFaulkned 
I1972D . Evidence for this picture comes from the metal- 
rich [Fe/H]~-h0.4 open cluster NGC 6791, where a 
t hird of its heliu m-burning stars are hot BHB stars 
( Yon g et al.l[200nh . 

It is unlikely, however, that a HVS is ejected during the 
horizontal branch phase. In the red giant phase, prior to 
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Fig. 4. — Distri bution of distance and velocity of simulated 1 
HVSs from the Br omle y et al. (2006) ejection model. The dashed 
line indicates the approximate escape velocity of the potential. In 
the range 10 < J? < 50 kpc, 95% of the simulated 1 Mq HVSs are 
bound, and 2/3 have crossed the Galactic center multiple times. 
If HVSs are BHB stars, we should thus find approximately equal 
numbers of stars at +300 km s~^ and —300 km s~^, contrary to 
the observations. 



becoming a BHB star, a star swells to a radius of ~1 AU. 
HVS ejections from a single M BH require stellar binaries 
with separations of 0.1 - 1 AU (|HillElll988trBromlev et aH 
|2006[) . Thus any progenitor HVS system experiences a 
common envelop e phase that in- spirals the binary system 
over < 10^ yrs (jWebbinkl |2007| ) that probably prevents 
a HVS ejection by the time the star reaches the BHB 
phase. If HVSs are BHB stars, they must be ejected as 
main sequence stars. This condition is relaxed for single 
stars ejected by a pair of MBHs or an in-spiraling IMBH. 

However, the lifetime constraint derived from the ob- 
served velocities means that the HVSs cannot be BHB 
stars descended from stars of mass less than 2 Mq. The 
progenitors of BHB stars span a wide range of stellar 
masses and include stars with main sequence lifetimes 
up to a Hubble time. Thus low-mass HVSs ejected many 
Gyr ago may be observed as BHB stars in our survey. In 
Figure S] we illustrate the spatial and velocity distribu- 
tion of such a popu lation from a simulation of 1 Mq HVS 
ejections. Bromlcv et al.l ([2006) describes the details of 
the ejection model. In brief, we disrupt equal mass bi- 
naries using the Hill's mechanism, and assume that the 
1 Mq stars are ejected at a random time during their 
main sequence lifetime. For reference, we also show the 
approximate escape velocity of the potential model with 
dashed- lines (Figure |4]). 

Our magnitude-limited survey detects BHB stars in 
the range 10 < i? < 50 kpc. Over this range of R, 95% 
of the simulated 1 Mq HVSs are bound; 2/3 of them 
have crossed the Galactic center multiple times. Thus, if 
HVSs arc BHB stars, we should observe approximately 
equal numbers of BHB stars at -1-300 km s~^ and —300 
km s~^, contrary to the observations (Figure [2|). "We 
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Fig. 5. — Location of our 7 HVSs (stars), the possibly bound 
HVSs (circles, triangles), and HVSl (plus sign) assuming they are 
main sequence B stars, z is the distance above the Galactic plane 
and r is the distance along the Galactic plane, such that R = 
(r^ -t-z^)""'^. For reference, we sketch the Milky Way and the Sun 
at (r, z) = (8, 0) kpc. 



simply don't see BHB stars from long-lived progenitors 
falling back onto the Galaxy with large negative radial 
velocities. 

Thus HVSs must have main sequence lifetimes less 
than the ~1 Gyr turn-around time. Given the color- 
selection of our survey, the HVSs are either 3-4 Mq 
main sequence stars, or possibly hot BHB stars evolved 
from M > 2 Mq progenitors. Our understanding of 
M > 2 Mq BHB stars is limited by the absence of stel- 
lar evolution tracks that connect the zero age main se- 
quence mass to the horizontal branch core mass. Main 
sequence stars, on the other hand, match the curr e nt un - 
derstanding of the Galactic center: iManess et al.l ()2007l ) 
find evidence for a top-heav y initial mass fun ction in the 
Galactic center, and Eisenhauer et al.l ()2005[ ) show that 
the stars presently in orbit around the central MBH are 
main sequence B stars, so me of which may be the forme r 
companions of the HVSs (iGinsburg &: Loebl[200l[200l . 
The 3-4 Mq main sequence stars sampled by our survey 
have lifetimes of 160 - 350 Myr, consistent with the ab- 
sence of bound HVSs falling back onto the Galaxy with 
large negative velocities. We conclude HVSs are likely 
main sequence stars. 

5. SPACE DENSITY OF HYPERVELOCITY STARS 

Our essentially complete survey places interesting con- 
straints on the space density of HVSs. The presence of a 
MBH in the Galactic center inevitably ejects a fountain 
of HVSs from the Galaxy. If HVS ejections are continu- 
ous and isotropic, the ir space density shou ld be propor- 
tional to (see also lBromlev et aLll2007[ ). The volume 
sampled by our survey is proportional to R^. Thus, in 
the simplest picture, we expect the number of HVSs in 
our survey to have the dependence iV(< R) oc R. The 
actual space density depends on the luminosities of the 
HVSs. We use the space density to predict HVS discov- 
eries for future surveys. 

5.1. HVS Luminosities and Distances 
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Fig. 6. — Observed distribution of HVSs (siars) vs. Galacto- 
centric radius K. Errorbars indicate 10% luminosity uncertainties. 
Dotted line shows the fraction of HVSs with travel times exceed- 
ing their lifetime at distance i?, indicated on the right-hand vertical 
axis. The best-fit line through the data [dashed line), corrected for 
incompleteness, has slope 0.159 it 0.017 kpc~^. 
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Fig. 7. — Cumulative number of 3 - 4 M0 HVSs vs. Galacto- 
centric radius R. Over 47r str, the space density of 3 - 4 Mq HVSs 
is p{R) = (0.077 ± 0.008)iJ-2 kpc'^. We infer there are 96 ± 10 3 
- 4 Mq HVSs within R < 100 kpc. 



Having established that HVSs are hkely main se- 
quence stars, we estim ate their intrinsic luminosities from 
iSchaller et al.l ()1992f ) stellar evolutionary tracks for solar 
metallicity stars. Five of the HVSs in our survey have col- 
ors (see Figure[T]) and spectra consistent with 3 Mq stars; 
the two bluest HVSs are consist ent with 4 Mp, stars. 
Using bolometric corrections from iKenvon fc Hartmannl 
(|1995h . we estimate My (3 Mq) ~ -0.3 and My (4 
Mq) ~ —0.9. These absolute magnitudes put the HVSs 
at distances in the range 45 < i? < 90 kpc. 

Figure[5]illustrates the HVS's location in the outer halo 
of the Galaxy. The possibly bound HVSs, in contrast, are 
mostly found at closer distances; there are 11 possibly 
bound HVSs and no unbound HVS in the bright FAST 
sample. Given a N{< R) cc R distribution, one may 
wonder why there are no unbound HVSs nearby. Pro- 
jection effects may provide one explanation. Consider 
SDSS J144955. 58+310351.4 (highlighted in Paper H), a 
possibly bound HVS located at i? ~ 17 kpc and traveling 
v,{ = +U7 km s-\ If SDSS J144955. 58-1-310351.4 has 
a purely radial trajectory, its true space velocity is -f 525 
km s~^ and it may thus be unbound. But because we 
do not know which of the possibly bound HVSs are true 
HVSs, we use only the well-defined sample of unbound 
HVSs below. 

We next compare the HVS distances with the range of 
distances sampled by our magnitude-limited survey. Our 
survey is complete over the magnitude range 15 < < 
19.5, corresponding to heliocentric distances of 12 < d < 
100 kpc and 16 < d < 130 kpc for 3 Mq and 4 Mg 
stars, respectively. The range of Galactocentric R varies 
with Galactic latitude and longitude on the sky. Over 
our 7300 deg^ survey region, we are complete over the 
range 20 < i? < 93 kpc and 24 < R < 124 kpc for 3 



Mq and 4 Mq stars, respectively. The inner edge of the 
MMT survey, in which the HVSs are found, is located at 
i? = 38 kpc and i? = 48 kpc for 3 Mq and 4 Mq stars, 
respectively. 

Figure [6] plots the cumulative distribution of distances 
for the 7 HVSs in this survey. The distribution is re- 
markably consistent with a continuous ejection process: 
a KS test finds a 0.96 likelihood of drawing the distri- 
bution of HVSs from a linear distribution A^(< R) oc R. 
This conclusion is independent of the absolute scale of 
R. 

5.2. Lifetime Correction 

Before computing the space density of HVSs from Fig- 
ure [6l we must correct for the number of HVSs missing 
in our survey bec ause of their finite life times. Using the 
ejection models of lBromlev et al.l ()2006[ ) . we calculate the 
fraction of stars with v > v^sc at a given R with travel 
times exceeding their lifetimes. The dotted line in Fig- 
ure [S] shows the weighted average of the 3 and 4 Mq 
correction factors, appropriate for our set of HVSs. 

The lifetime correction is negligible for 3 Mq HVSs 
in our magnitude range, but it is substantial for the 
faintest 4 Mq HVSs. The two 4 Mq HVSs in our sur- 
vey are located at i? = 56 kpc and i? = 82 kpc; given a 
N{< R) (X R dependence, we expect to find one or two 
more 4 Mq HVSs ou t to i? = 124 kpc. However, the 
iBromlev et al.l (|2006f ) ejection models show that 50% of 
unbound 4 Mq HVSs are dead by i? = 90 kpc, and 67% 
are dead by i? = 110 kpc. Thus the absence of 4 Mq 
HVSs with i? > 90 kpc is consistent with the interpre- 
tation that they are main sequence stars with short, 160 
Myr lifetimes. The 3 Mq HVSs, on the other hand, fiU 
the region where we can detect them, consistent with 
their longer, 350 Myr main sequence lifetimes. 
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5.3. Density Distribution 

We derive the HVS space density from a least-squares 
fit to the observations, corrected for lifetime. The slope is 
0.159 ±0.017 kpc~^ (see Figure^]). We estimate the un- 
certainty by assigning random uncertainties to the HVSs, 
drawn from a Gaussian distribution with a =10% in lu- 
minosity, and by re-fitting the line 10* times. Because all 
7 HVSs have g'g > 17, they are drawn from the MMT sur- 
vey and cover an effective area of 6800 deg^ of sky (93% 
of 7300 deg^). Corrected to in steradians, the space den- 
sity of 3 - 4 Mo HVSs is thus p{R) = (0.077 ±0.008)i?-2 
kpc~^. 

Figure[7]plots the cumulative number of 3 - 4 Mq HVSs 
in the Galaxy. Our observations imply there are 96 ± 10 
late B-type HVSs in the sphere R < 100 kpc, about 1 
HVS per kpc. 

5.4. Prediction for Future Surveys 

We can use the space density of 3 - 4 Mq HVSs to 
predict HVS detections in future surveys. Our predic- 
tions make the following assumptions. (1) We assume 
that the distribution of HVS velocities is the same for 
all stellar masses. Ejection velocity depends weakly on 
the mass of the stellar binary (mi -I- m2y^^ in the Hill's 
mechanism, but there is no such dependence for stars 
ejected from a binary MBH. (2) We use the Salpeter 
mass function, M~'^-^^ (|Salpeteij|1955f ). and the present- 
day mass function of the b ulge, M~ for M > 1 Mq 
and M-2-35 f^j. m < 1 M© (jMezger et al.lll99 9'). as two 
representative mass functions for calculating the number 
ratios of stars of different stellar masses. (3) If HVSs are 
ejected by the Hill's mechanism, we implicitly assume 
that the same fraction of stars are in compact binaries 
at all masses; no such assumption is re quired for the 
binary MBH mechanism. (4) We use the iGirardi et al] 
(|200l stellar isochrones of solar abundance to obtain lu- 
minosities for stars at a given mass, or color. We use 
the luminosities to calculate survey volumes for a set of 
magnitude limits. (5) Finally, we assume observation of 
100 % of stars in a giv e n area of sky. 

iKollmeier fc Gouldl (pOOTh propose observing faint 
19.5 < g' < 21.5 stars near the main-sequence turn-off 
0.3 < {g' — i') < 1.1 as the optimal strategy for finding 
low-mass HVSs. The {g' -i') colors select 0.8 < M < 1.3 
Mq solar metallicity stars with absolute magnitudes 
3.7 < My < 6.3. Using the present-day mass func- 
tion of the bulge, we estimate that there is 1 HVS per 
~50 deg^ in the proposed magnitude range, in excel- 
le nt ag reemen t with the 1 H VS per 45 deg^ estimated 
bv IKollmeier fc Gouldl ()2007l ). However, a Salpeter mass 
function predicts an order-of-magnitude lower density, 
1 HVS per ~500 deg^. The ratio of high- to low-mass 
HVSs thus provides a sensitive me asure of the stellar 
mass function near t he cen tral MBH (jKollmeier k Gouldl 
HOOT; Bro mlev et aI1l200l . 

SEGUE is an on-going survey that includes spec- 
troscopy of 114 4 stars along ■-■^200 sightlines w ith the 
SDSS telescope (jAdelman-McCarthv et alll2007bl ). Each 
sighthne covers 7 deg^ of sky, for a total of 1400 deg^ 
of spectroscopic coverage. Notably, SEGUE acquires 
spectra for 150 BHB/A-type stars per sightline over 
the magni tude range 14 < g' < 20.5. T he BHB/A 
color cut (jAdelman-McCarthv et all l2007"bl ) selects for 




120 



Fig. 8. — Minimum rest-frame velocity and distance of the HVSs 
(stars) and HVSl (plus sign). Errorbars indicate the precision of 
the measurements; systematic errors may be larger. Travel times 
from the Galact ic center (dot ted lines) are calculated using the po- 
tential model of lBromlev et al. (2006) and assuming the minimum 
rest frame velocity v^f is the full space motion of the stars. The 
distribution of travel times statistically favor a constant ejection 
rate, though a 120 Myr-old burst of HVSs is also allowed. 



1.5 < M < 2.5 Mq stars with absolute magnitudes 
0.7 < My < 2.8. We estimate SEGUE should find -5 to 
~15 HVSs for the Salpeter and present-day mass func- 
tions, respectively. However, if only a fraction of BHB/A 
stars are observed over the 1400 deg^ area, or if the spec- 
tra are not reliable to a depth g' — 20.5, our predictions 
may be over-estimates. 

SEGUE also targets F/G-type stars over the magni- 
tude range 14 < < 20.2. If we ignore low- metallicity 
targets as unlikely HVS candidates, SEGUE targets 50 
F/G main sequence s tars per sightline. The color cut 
0.2 < {g' - r') < 0.48 (lAdelman-McCarthv et al']l2007bD 
selects for 1.1 < M < 1.4 Mq stars with absolute mag- 
nitudes 3.0 < My < 4.6. We estimate that 50 stars 
per sightline account for ^10% of the total population of 
F/G stars, based on number counts of stars in stripe 82. 
We estimate that SEGUE should find ~0.1 to -1 F/G- 
type HVS for the Salpeter and present-day mass func- 
tions, respectively. The F/G-type HVS discovery rate is 
low because SEGUE is not complete for F/G stars and 
because F/G stars have intrinsically low luminosity. 

Regardless of stellar type, future HVS discoveries re- 
quire surveys that 1) sample large areas of sky » 10^ 
deg^ and 2) sample large depths i? » 10 kpc. As a prac- 
tical matter, the HVSs must also have a high contrast 
with respect to the indigenous stellar population for ef- 
ficient detection. Our HVS survey is successful because 
late B-type stars are luminous, and because the contam- 
ination from evolved BHB stars and their low mass pro- 
genitors is minimal, ^100 halo stars per HVS. 

6. EJECTION HISTORY OF HYPERVELOCITY STARS 
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Knowing the distances and velocities of the HVSs, we 
can estimate their travel times from the Galactic cen- 
ter. Although we do not know the HVSs' true space mo- 
tions, at large distances their motion should be mostly ra- 
dial. Thus we calculate travel times assuming the HVSs' 
minimum rest-frame velocities are their full space mo- 
tions. This assumption means the travel time estimates 
ar e in fact upper limits . We use the potential model 
of iBromlev et al.l (|2006f ) to correct for the deceleration 
of the HVSs. The deceleration correction reduces travel 
times by 10% compared to travel times calculated assum- 
ing constant velocity (as reported in Paper I) . Table [1] 
Hsts travel time estimates for HVS4 - HVSIO and HVSl, 
though we do not include HVSl in this analysis. Figure 
[8] plots the distribution of travel times. 

The HVSs sample the ejection history of stars from 
the Galactic center over the past '^165 Myr. The HVS 
travel times might favor a single burst, like that expected 
from the in-spiral of a binary MBH, or a continuous ejec- 
tion process, as expected from Hill's mechanism. Inter- 
estingly, 5 of the 7 HVS travel times appear clustered 
around 120 Myr. A Kolmogorov-Smirnov test finds a 0.8 
likelihood of drawing all 7 HVSs from a linear distribu- 
tion in time, and a 0.2 likelihood of drawing them from 
a single burst (e.g. travel times that fall along the dotted 
line marked 120 Myr in Figure [8|). 

The statistics appear to favor a continuous ejection 
process. Thus Hill's mechanism can plausibly account 
for the entire set of observed HVSs. A single burst event 
is also allowed, however. If HVSs were ejected from the 
in-spiral of a binary MBH, the observations suggest such 
an event may have occurred ~120 Myr ago. We caution 
that small number statistics prevent any strong conclu- 
sion; additional HVS discoveries are required to properly 
constrain the ejection history of stars from the Galactic 
center. 

7. CONCLUSIONS 

Our targeted HVS survey, a spectroscopic survey of 
stars with B-type colors, is 96.5% complete over 7300 
deg^ of the SDSS DR5. We report 3 new HVS discoveries, 
for a total of 7 HVSs with 15 < < 19.5 over a sixth 
of the sky. In addition, we find 19 possibly bound HVSs 
with velocities 275 < Wrf < 400 km s^^. 

The significant excess of +300 km s^^ stars compared 
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to —300 km s^^ stars shows that the HVSs must be short- 
lived. Any population of post-main sequence HVSs must 
contain stars falling back onto the Galaxy, contrary to 
the observations. We conclude the HVSs are likely 3-4 
Mq main sequence stars, although hot BHB stars from 
M > 2 Mq progenitors are also a possibility. Observa- 
tions of stars in the central 10 pc of the Galaxy are now 
probing ^3 Mq main sequence stars, some of which may 
feed the MBH and produce HVSs. Thus there is a need 
to better understand the BHB phase of 2 - 3 Mq stars 
in the context of HVSs and Galactic center research. 

The spatial distribution of HVSs supports the main 
sequence interpretation. The longer-lived 3 Mq HVSs 
fill our survey volume; the shorter-lived 4 Mq HVSs are 
missing at faint magnitudes. The spatial distribution 
of HVSs is remarkably consistent with a A^(< R) oc R 
distribution, as expected for a continuous ejection pro- 
cess. The density of 3 - 4 Mq HVSs is p{R) = (0.077 ± 
0.008)i?-2 kpc-3, implying a total of 96 ± 10 such HVSs 
within R < 100 kpc. 

A sample of > 100 HVSs is required to discriminate 
among HVS ejection mechan isms and Galactic potential 
models (jSesana et al.ll2007a^ . Thus enough 3 - 4 Mq 
HVSs may exist in the Galaxy to test the ejection and 
potential models, especially if the bound HVSs can be 
used. SEGUE and other future surveys will find more 
HVSs and will allow further use of this important class 
of objects to connect the Galactic center to the Galactic 
halo. 
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TABLE A2 
HVS SURVEY: NEW OBSERVATIONS 



RA 


Dec 


9' 


{«' - 9')o 




WD 


Vq 


J2000 


J2000 


mag 


mag 


mag 




km s ^ 


0:23:09.05 


-0:33:42.0 


16.277 ±0.020 


0.250 ±0.028 


-0.306 ± 0.039 


1 


-13 ± 31 


0:23:53.29 


-1:04:46.4 


18.304 ±0.016 


0.749 ± 0.042 


-0.255 ±0.025 





20± 10 


0:43:50.55 


-9:52:27.0 


16.326 ±0.032 


0.832 ±0.038 


-0.358 ±0.038 





-117 ± 11 



Note. — Table K2l is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here 
for guidance and content. 
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APPENDIX 
DATA TABLE 

Table IA2I [available as tab2.tex in the source] lists the 381 new observations from our HVS survey, excluding the 
extra-galactic objects. 206 HVS candidates were observed with the MMT and 175 were observed with FAST (see §2). 
Table [X2] includes columns of RA and Dec coordinates (J2000), g' apparent magnitude, [u' — g')o and {g' — r')o color, 
and our heliocentric velocity vq. The column WD indicates whether the object is a white dwarf (WD=1) or a B-type 
star (WD=0). 



